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Abstract. We describe the construction and operation of a time-orbiting potential 
trap (TOP trap) that has different oscillation frequencies along its three principal axes. 
These axes can be rotated and we have observed Bose-Einstein condensates of ^^Rb 
with a rotating ellipsoidal shape. Under these conditions it has been predicted that 
quantized vortices form and are stable. 
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1. Introduction 

The first experimental observations of Bose-Einstein condensation in alkali metals in 
1995 |]T], 0] and subsequent experiments have led to much progress, both theorerical 
and experimental, in understanding the physics of these condensed systems (for a 
review see f^). Initial experiments focussed on the macroscopic properties of this new 
state of matter, such as collective excitations |^. Later experiments on atom lasers 
P, and the production of mult i- component condensates ^ have focussed on the 
coherence properties of EEC's. In particular the description of a condensate by a single 
wavefunction with a definite phase has been beautifully demonstrated by matter-wave 
interference experiments 0. Nevertheless, in contrast to earlier work on liquid 



helium |jT^ , the spectacular phenomena of superfluidity have yet to be comprehensively 
demonstrated and studied in the alkali metal condensates e.g. non-viscous flow in narrow 
channels, quantised vortices and persistent currents. Some of these phenomena are a 
consequence of the irrotational nature of superfluid flow, which in turn stems from 
the existence of the single-valued wavefunction. Recently much theoretical effort, by 



various groups [|I2], |I3|, |T^, 0, has been devoted to possible methods of formation 
and observation of quantised vortices and persistent currents, which are associated with 
rotation of the condensate. It has been shown that rotating an anisotropic trap should 



create vortices in a similar manner to those observed in liquid helium [|1^ . 



In contrast to the time-averaged orbiting potential (TOP) trap used in this work, 
it is difficult to create a rotating potential using a magnetostatic trap such as the 



loffe-Prit chard |]T8l configuration. However it may be possible to study superfiuidity 
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using other methods that have recently been suggested: stirring a trap with hght forces 
||T9| , imprinting a 27r phase-winding with hght beams and manufacturing a vortical 



wavefunction in a spinor superposition state [21 1. Only the latter method has been 



successful so far; the challenge of observing a vortex for a condensate in a single internal 
state has yet to be met. 

However these alternative methods cannot directly test some of the interesting 
theoretical predictions which have recently been made for states in a rotating trap which 
have non-integer values of angular momentum (in units of Plancks constant) per particle 

The method presented here allows rotation without changing the ellipticity of the 



potential even for large ellipticities. These features cannot be achieved in the modulation 
method [HI used to excite m = 2 condensate excitations. 



2. Experimental arrangement and BEC in TOP trap 

We use a combination of laser cooling, magnetic trapping and evaporative cooling to 
obtain a ^''Rb condensate in the F = 2,mF = 2 state in a similar manner to the first 
observation of BEC U|. The experimental arrangement is similar to those described 



elsewhere p3| , |2^ , and only important features of our setup are given below. 
The light for our experiment comes from two external cavity diode lasers, master and 
repumper, operating close to the frequency of the rubidium D2-transitions. The light 
from the master laser is injected and amplified by a tapered semiconductor amplifier 
chip in a master oscillator power amplifier (MOPA) |^ configuration which provides 



~ 500 mW of output. The light is fed to two magneto-optical traps (MOT) in a double 
trap configuration [^] to preserve a good vacuum in the experimental region while 
obtaining high loading rates. 

A novel feature is the use of a pyramidal configuration of mirrors to form the 



collection MOT ||2^. A large number of rubidium atoms are emitted from a getter 
source [28| in this region and then continuously cooled and transferred down a 30cm 
pipe of 15mm inner diameter, into the region of higher vacuum where they are captured 
in the second MOT. Magnetic guiding, in the transfer pipe, increases the flux of atoms 
to the second MOT by a factor of 2 [^]. The second MOT is a standard 6-beam MOT 
configuration used to capture and store the atoms. In order to obtain large numbers of 
atoms we found that it was beneficial to keep the small (~ 2G) time-varying magnetic 
field used in the TOP trap switched on during loading into the second MOT. The 
Zeeman shift of this field causes atoms near the centre to be detuned from resonance, 



giving a partial dark MOT . 

After collecting 5 x 10^ atoms in the second MOT they are compressed, cooled to 
a temperature of 30 in an optical molasses and then optically pumped into the 
F = 2,mF = 2 state. The quadrupole field is created by two coils in anti-Helmholtz 
configuration with 500 turns each, operating at a maximum current of lOA. The 7kHz 
rotating field is created by two sets of Helmholtz coils with 5 turns in each coil. The 
current to each set of Helmholtz coils is supplied via matching transformers from a 
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600Vr audio amplifier. 

Once the atoms are loaded into a mode-matched TOP trap {25G bias field and 
65G/cm radial gradient) we use the following procedure to obtain BEG (see also Fig.|l| 
(a)). Firstly the TOP and quadrupole fields are adiabatically ramped to 40G and 
lOOG/cm in 2s to increase the collision rate. In a combined evaporation and compression 
ramp the radial quadrupole field gradient continues to rise to 194G/cm while the TOP 
field decreases to 20G in 20s. Now the TOP field itself is used to evaporate the hottest 
atoms by cutting to IG in 28s obtaining a final radial trap frequency of 175if2||[ Finally 
a 5s radio frequency cut, finishing 60% of the way in from the radius of the locus of 
-8 = {i' final = 0.98MHz), is used to obtain a pure condensate of ~ 3 x 10'' atoms. In 
Fig.|l](a) it is possible to follow our path to BEG for ^^Rb in a phase space diagram and 
compare it with the case of ^^^Gs that has also been studied in a parallel experiment 



[using the same magnetic trap) at our lab ||31 



The pictures in Fig.|l|(b) show the emergence of the pure condensate as the depth of 
the final RF cut is increased. We observe the sudden increase of optical depth and the 
reversal of shape characteristic for a BEG in time of flight (TOF) imaging. A laser 
pulse of 10/is, synchronised with the TOP field and resonant with the F = 2, nip = 
2 — s> F' = 3, m'p = 3 transition, is used for absoption imaging of the cloud of atoms in 
the xz-plane. 

3. Elliptical TOP potential 



The invention of the TOP trap |2^ lead to the first observation of BEG |1|] in a dilute 
alkali vapour of rubidium, and this type of trap has subsequently also been used by other 



groups to trap caesium ^ and sodium IQ. The TOP trap is formed by combining 



a spherical quadrupole field, which can be written as 

Bq = B'q {xe^ + yey - 2ze^) (1) 
near the origin [B'q is the gradient), with a rotating field in the xy-plane 

B = 5^ cos Uot + By sin ujQt (2) 

The rotation frequency Uq is chosen to be greater than the oscillation frequencies of 
atoms in the trap Ux, ojy and uj^., but less than the Larmor frequency so that atoms do not 
change m^-state. The TOP trap is usually operated in the axisymmetric configuration 
where the amplitudes of the fields in the two directions are equal, B^ = By, and the 
trapping frequencies are in the ratios cux = cOy = tOz/ ||. However it is numerically 
straightforward to calculate the potential in the more general case where the ratio of the 
magnitudes of the rotating fields is B^/By = E. The TOP potential is the time-average 
of 

I The trapping frequencies were measured by exciting 'dipole' oscillations of a small thermal cloud 
with a sudden change of the bias magnetic field. 

§ A non axisymmetric TOP trap has been used to obtain BEG in Na M. 
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U{x, y, z, t) = /x-BqKx + Evq cosco'ot)ea; + {y + tq smuQt)ey — 2zez\ (3) 

where /i is the magnetic moment of the atoms and vq = By/ B'q is the distance from the 
centre of rotation at which the locus 5 = intersects the y-axis. For E = 1 this locus 
is a circle of radius tq in the xy-plane, and in general it is an ellipse of eccentricity E in 
the xy plane described by the position vector 

r = To (-E cos ujot e^. + sin uj^t By) (4) 

Without loss of generality we can consider ii^ > 1 so that the order of the numerically 
calculated trapping frequencies, shown in Fig.^, is uj^ < oJy < uj^. The harmonic 
approximation to the TOP potential 

U = Un(ulx^ + uy + ujlz^) (5) 

is accurate at distances small compared to tq. We are mainly interested in the 
deformation of the time-averaged potential in the xy-plane, characterised by the 
parameter Uy/ux = e > 1. Contours of constant energy in the x?/-plane are ellipses 
given by 

+ e^y'^ =constant (6) 

whose major axes are in the same direction as the major axis of the ellipse in Eq. 
Note however that the eccentricity of the TOP is much less than E and that the field 
in the x-direction reduces all the trap frequencies, not just Ux- Thus it is not possible 
to reach very large deformations of the TOP potential. For small deformations the 
ellipticity in the trap frequency varies as de/dE = 1/4. 



4. Rotation of the TOP 



A rotation of an angle (p is mathematically obtained by a transformation of coordinates 
expressed by the matrix 

cos (f) sin (f) 
— sin (j) cos (f) 

Applying this with a time varying angle (p = VLt to Eq.^ gives the position vector of a 
rotating ellipse 

r' = To (-Ecosfit coscuo^ + sini7t sincjot) 63;+ , , 

To (— i^sin Vtt cos ujQt + cos fit sint^ot) e^. 

When E = 1^ this expression reduces to sin(co'o — fi)t and cos(co'o — fi)t, which corresponds 
to the usual TOP configuration with a very slightly different frequency, since f2 <^ cuq. 
(e.g. fi < 27r X lOOifz whereas c^o = 27r x 7kHz in these experiments). To make the 
locus of -B = a rotating elhpse as described by Eq.|, we need to apply magnetic fields 
in the x and y directions which vary in time in the way described in these equations. 
This is achieved using the circuit shown schematically in Fig.|^, starting with four signals 
proportional to Ecosuot, sincuot, cos fit and sin fit. These are then multiplied in pair 
by four voltage multiplier to give the four product voltages. 
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A = cos fit E cos (Jot 

C = sin fit sinuot , . 

B = sin nt E cos cuot 
D = cos fit sincjo^ 

These voltages are then added (or subtracted) by simple operational amplifier circuits 
to yield final signals 

oc A + C 

By oc D-B, ^ ^ 

which are amplified by the audio amplifier and fed to the coils producing the time- 
varying fields along x and y. 

It is important to ensure that any high frequency noise introduced in the multiplication 
and summing stages is filtered out before the amplifier otherwise it leads to a large 
decrease in the lifetime of the trapped atoms. This 'noise-induced' loss becomes more 
severe at lower values of the rotating bias field, and our present system cannot be 
operated below IG. 



5. BEC in the rotating anisotropic TOP trap 

A simple description of the behaviour of trapped atoms in the elliptical rotating trap, 
is obtained from the classical equations of motion for anisotropic harmonic oscillation 
of a particle in a frame rotating at angular frequency fi about the z-axis 

X = - (cj2 - fi2) a; - 2fi y 

y = -{u;'y-n^)y + 2nx (11) 

z = z 

These equations are discussed in many mechanics textbooks e.g. in the book by Lamb 
||33|| , he shows that the two equations for the xy motion describe Blackburn's pendulum, 
a pendulum with different effective lengths in two orthogonal directions on a turntable. 
The oscillation frequencies in the x and y directions are reduced by the effect of the 
centripetal acceleration outwards and the x and y equations are coupled by the Coriolis 
force f2 X f . One obvious feature of these equations is that the motion becomes unstable 
when the rotation is equal to the lowest trap oscillation frequency, fi = u;^;. 

We use the procedure outlined in section 2 to obtain a BEC in a normal circular 
TOP trap and then adiabatically change the amplitude of B^-, By or both to the 
desired ellipticity E. We found that the BEC remained condensed during this adiabatic 
expansion, during which the RF knife was switched off. Fig.^ shows time of flight 
images of the condensate after release from the rotating elliptical trap. By releasing 
the BEC when the probe beam was aligned with its long or its short axis the evolution 
of the X to z and y to z aspect ratios of the BEC was monitored. The inversion of 
the shape as the time of flight increases clearly demonstrates that the expansion of the 
cloud is dominated by its self energy. For the data shown in Fig.^ B^ remained at its 



Bose-Einstein condensation in a rotating anisotropic TOP trap 



6 



final vaue of 2G and By was ramped up to a value of 8G in 2s giving an ellipticity of 
the time-averaged potential of E = 4 and e ~ 1.4. These fields imply a static trap 
frequency of 70Hz. With a rotation frequency of 30Hz, one obtains an effective trap 
frequency (cu^ — = 63Hz. We numerically solved the Castin-Dum |]3^ equations 

for the effective trap frequencies in a non-rotating frame. Fig.|^ shows the comparison of 
experimental and theoretical data, which are in good agreement. This proves that our 
BEC has survived the rotation of the elliptical trap and that the Coriolis effects have 
had negligible influence. 



6. Conclusions 



We have constructed an anisotropic and rotating TOP trap and shown that a BEC may 
be created and transferred to the rotating trap. The measured properties at low rotation 
frequencies agree well with the existing theory. Faster rotation will be a useful tool for 
further studies of superfiuid properties of a BEC, such as nucleation and stabilization 
of vortices. 
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Figure 1. (a) Plot of phase space density vs. number of atoms during the evaporative 
cooling ramp. We compare the evaporative cooling efficiency for *^Rb and and ^^'^Cs 
in two identical magnetic traps. "Crunch" refers to the combined evaporation and 
compression ramp, (b) BEC formation in a trap with IG bias field. The images were 
taken after 7ms TOF using absorption imaging. 
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Figure 2. Plot of the trapping frequencies in an elliptical TOP trap as a function 
of the ellipticity E = B^/ By where B^ is varying with E for fixed values of By and 
the quadrupole gradient. The frequencies are plotted as fractions of their values for 
Bx = By. The dotted line (uz) remains approximately equal to the mean of the other 
frequencies {lOx +'^j/)/2. 
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Figure 3. Schematic of the circuit for modulating the sine and cosine waves driving 
the TOP coils. 
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Figure 4. Pictures of BECs at different times of flight after being released from the 
elhptical rotating TOP trap. The top row shows BECs released with their x axes 
parallel to the imaging axes. The bottom row shows BECs released with their y axes 
parallel to the imaging axes. 



Bose-Einstein condensation in a rotating anisotropic TOP trap 



12 



4.0 




' 1 ' 1 ' T" 

5 10 15 

Time [ms] 




Figure 5. The expanison of the BEC is shown as a function of time. The aspect ratios 
evolve under the self-energy expansion. The theoretical curves have been obtained 
solving the Castin-Dum equations without free parameters, (a) Evolution of the x/z 
aspect ratio, (b) Evolution of the y/z aspect ratio. The data points are the average of 
four shots and the error bars represent their shot to shot variation. 



